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Abstract 

In plant leaves, resource use follows a trade-off between rapid resource capture 
and conservative storage. This "worldwide leaf economics spectrum" consists of 
a suite of intercorrelated leaf traits, among which leaf mass per area, LMA, is 
one of the most fundamental as it indicates the cost of leaf construction and 
light-interception borne by plants. We conducted a broad-scale analysis of the 
evolutionary history of LMA across a large dataset of 5401 vascular plant spe- 
cies. The phylogenetic signal in LMA displayed low but significant conserva- 
tism, that is, leaf economics tended to be more similar among close relatives 
than expected by chance alone. Models of trait evolution indicated that LMA 
evolved under weak stabilizing selection. Moreover, results suggest that different 
optimal phenotypes evolved among large clades within which extremes tended 
to be selected against. Conservatism in LMA was strongly related to growth 
form, as were selection intensity and phenotypic evolutionary rates: woody 
plants showed higher conservatism in relation to stronger stabilizing selection 
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and lower evolutionary rates compared to herbaceous taxa. The evolutionary 
history of LMA thus paints different evolutionary trajectories of vascular plant 
species across clades, revealing the coordination of leaf trait evolution with 
growth forms in response to varying selection regimes. 



doi: 10.1002/ece3.1087 



Introduction 

Although leaf morphology has evolved manifold varia- 
tions across vascular plant species, there is strong evi- 
dence of a universal spectrum constraining leaf 
functioning from rapid resource capture to efficient 
resource use (Grime et al. 1997; Reich et al. 1997; Diaz 
et al. 2004; Wright et al. 2004). This worldwide leaf eco- 
nomics spectrum (Wright et al. 2004) consists of the cor- 
related variation among several key plant traits, including 
leaf mass per area (LMA, the ratio between leaf dry mass 
and leaf area) which describes the dry mass investment 
for light interception per unit leaf area (Lambers and 
Poorter, 1992). LMA captures a central axis of functional 
variation in plants (Grime et al. 1997; Grime 2001; West- 
oby et al. 2002; Diaz et al. 2004), correlating negatively 
with mass-based photosynthetic rate and leaf macronutri- 
ent concentrations (Reich et al. 1997; Wright et al. 2004), 
and positively with leaf life span (Wright et al. 2004). At 
large spatial scales, LMA varies across climatic gradients, 
displaying on average higher values in hotter, drier, and 
higher irradiance habitats, particularly once other factors, 
such as deciduousness or plant functional type composi- 
tion, are taken into account (Reich et al. 1997; Wright 
et al. 2004; Reich et al. 2007; Poorter et al. 2009). 

Despite a sound knowledge of the physiological and 
ecological correlates of LMA, we remain largely ignorant 
of the evolutionary history that gave rise to the present- 
day wide variation in this key trait. Two types of overall 
evolutionary behavior may be expected for a vegetative 
trait such as LMA: it can either be evolutionary labile and 
vary independently from the phylogeny across species, 
which can lead to high functional convergence across 
clades (Kraft et al. 2007), or it may alternatively display 
phylogenetic patterns structured by selection within plant 
lineages (Blomberg et al. 2003; Losos 2008). In the labile 
trait hypothesis, we nonetheless expect consistent evolu- 
tion with other plant traits, or syndromes such as growth 
form. Recent advances have indeed suggested that the 
evolution of growth form, as a syndrome of multiple 
traits, has been a major driver of functional trait evolu- 
tion (Moles et al. 2005; Kerkhoff et al. 2006). Marked 
taxonomic patterns in LMA would therefore relate more 
strongly with growth forms than with phylogenetic relat- 
edness. Alternatively, following the second hypothesis. 



LMA would have evolved under strong selection and be 
conserved along plant lineages, which would generate a 
consistent and detectable phylogenetic signal. 

Phylogenetic signal can be broadly defined as the infor- 
mation conveyed by the variation in phenotypic trait val- 
ues within and across clades along a phylogeny. In the 
present study, we conducted a comprehensive analysis of 
the evolutionary history of leaf economics across the phy- 
logeny of vascular plants based on the analysis of the phy- 
logenetic signal in LMA. Variation in trait similarity 
among taxa may reveal patterns of either trait diversifica- 
tion or trait conservatism within clades, that is, the ten- 
dency for relatives to be more similar than expected by 
chance alone. Such macroevolutionary patterns may pro- 
vide insight into the microevolutionary processes (e.g., 
selection and drift) shaping trait evolution within lineages 
(Hansen and Martins, 1996; Hansen 1997; Diniz-Filho 
2001). 

Variation in leaf mass per area was charted across a dated 
phylogeny of 5401 species spanning all major clades of vas- 
cular plants. Using phylogenetic comparative analyses, we 
analyzed the phylogenetic signal in LMA and compared 
observed patterns to expectations based on different models 
of trait evolution. We show that LMA exhibits low but sig- 
nificant overall phylogenetic conservatism across vascular 
plant clades. We identify different clades displaying signifi- 
cant phylogenetic patterns of either trait conservatism or 
diversification. Moreover, we tested the interaction between 
the evolutionary histories of LMA and growth form across 
species. We note that potential additional explanatory vari- 
ables such as biogeography or pedo-climatic conditions 
were outside the scope of the present study. We found that 
patterns of LMA evolution consistently differed across 
growth forms, revealing evidence for a higher conservatism 
and slower trait evolution in woody species than in herba- 
ceous species. 

Material and Methods 
Trait data 

We collected LMA data for vascular plants {Tracheophyta) 
from 180 published and unpublished studies and elec- 
tronic databases (see Appendix SI). We only retained data 
collected in the field for outer-canopy leaves measured 
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following standardized protocols (Cornelissen et al. 2003) 
and for adult plants, to limit ontogenic effects. Mean val- 
ues were calculated for species with multiple records. The 
sample set of species and sites represented a wide range 
of plant communities in most of the climates where vas- 
cular plants occur, from Arctic tundra to tropical forest, 
from hot to cold deserts, and from grassland to wood- 
land. In total, we obtained LMA values for 5401 species 
in 241 families and 1835 genera: 5239 Angiospermae, 81 
Gymnospermae, 74 Monilophyta, and 7 Lycopodiophyta 
(see Fig. 2 for an overview of sample completeness). We 
checked that the variability in LMA occurred mostly 
among species compared to within species (across sites): 
^among = 16286 vcrsus Vwithin = 1912. In order to limit 
the influence of extreme values, we performed compara- 
tive analyses on log-transformed LMA values that we 
hereafter note I. 

Dated phylogenetic supertree 

Phylogenetic relationships between species were described 
as an informal supertree based on published phylogenies 
(Bininda-Emonds 2004). We used the tree of Chaw et al. 
(2000) as a backbone for relationships between major vas- 
cular plant clades and between families within Gymno- 
spermae on which we branched family-level trees for the 
Angiospermae (Davies et al. 2004), Monilophyta and 
Lycopodiophyta (Wikstrom and Kenrick 2001; Smith et al. 
2006). We delimited major clades Tracheophyta accord- 
ing to the phylogenetic nomenclature of Cantino et al. 
(2007). Family names in Angiospermae were matched to 
the latest phylogeny by the Angiosperm Phylogeny Group 
(APG III 2009). Within 34 large Angiospermae families, 
we resolved relationships among genera using published 
family phylogenies (see Appendix SI). Species were 
branched as polytomies within genera. The resulting 
supertree had 1675 internal nodes among which 58% were 
resolved as bifurcations (nodes with two daughter clades). 
Below the genus level, the percentage of bifurcations 
among ancestral nodes rose to 75%. In order to estimate 
branch lengths, we dated the ancestral nodes of the super- 
tree in a two-step procedure. We first dated well-identified 
nodes in the supertree using published ages. A total of 187 
ancestral nodes matched dated nodes from the literature 
(11% of internal nodes). For clades in Angiospermae, ages 
were taken from a comprehensive update of divergence 
times based on the analysis of sequence data (Bell et al. 
2010). We completed these ages with estimates derived 
from several sources for non-angiosperm clades (Wikstrom 
et al. 2001; Bremer et al 2004; Janssen and Bremer, 2004; 
Anderson et al. 2005). Node ages ranged from 4.4 Myr for 
the Juglandaceae family to 535 Myr for the tree root, 
namely the Tracheophyta divergence (mean: 88.1 Myr SD: 



75.3 Myr). Second, node ages were estimated using the 
bladj module of the Phylocom software (Webb et al. 2007). 
This procedure sets branch lengths by placing the nodes 
evenly between dated nodes, and between dated nodes and 
tips (of age 0). This has the effect of minimizing variance in 
branch length, within the constraints of dated nodes (Webb 
et al. 2007). The phylogenetic distance between two species 
was estimated by the age of their most recent common 
ancestor. 

An alternative method would have been to reconstruct 
a phylogeny from molecular data. For large datasets, how- 
ever, it is difficult if not impossible to obtain resolved 
phylogenies at species level. A recent study compared 
phylogenetic patterns in phenology for ca. 4000 species 
(Davies et al. 2013) using phylogenetic analyses of both 
an informal supertree, such as the one constructed here, 
and a phylogenetic tree obtained from molecular data 
and resolved at genus level. The patterns described were 
similar and led to the same conclusions. We also explored 
this option, but were only able to locate genetic data for 
about 40% of species. We eventually chose to retain a 
much wider coverage but lower resolution. 

Models of continuous trait evolution 

We tested and compared alternative evolutionary hypoth- 
eses by fitting different evolutionary models to our data. 
Models of trait evolution describe the macroevolutionary 
patterns which would be expected from hypothetical 
microevolutionary processes (Hansen and Martins, 1996). 
First, we considered a simple model of Brownian motion 
evolution (BM), which assumes that the trait evolves 
independently in each lineage by means of genetic drift 
and/or directional selection under environmental condi- 
tions fluctuating randomly and rapidly compared to evo- 
lutionary time (Hansen and Martins, 1996; Diniz-Filho 
2001). The BM model supposes evolution by random 
motion of trait values at a constant rate a along the 
branches of the phylogenetic tree and thus mimics the 
effects of drift (Hansen and Martins, 1996; Freckleton 
et al. 2002): during an infinitesimal period dt, the varia- 
tion in the trait value I is dl = a^dt, where parameter a 
controls the magnitude of stochastic perturbations during 
the course of evolution, or drift (Hansen 1997). 

Second, we considered models of LMA evolution incor- 
porating both stabilizing selection and drift (Hansen 
1997). These models conform to evolution following an 
Orstein-Ulhenbeck (OU) process: during an infinitesimal 
period dt, the variation of the trait value sums as the 
effects of drift, a^dt as in BM, and selection toward a 
phenotypic optimum 0: dl = —a.{l—6) + a~dt, where the 
parameter a controls the rate of adaptive evolution to 
the optimum (Hansen 1997). Parameter a measures the 
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magnitude of a supposed selective force. Note that a null 
value of a. leads to the BM model, which thus appears as 
a special case of OU model. Interestingly, OU models 
allow for different selective optima to be specified within 
different clades along a single phylogenetic tree (Hansen 
1997; Butler and King, 2004): changes in the selective 
optimum mimic variation in the selection regime along 
the phylogeny. After a change in the selection regime, 
evolution unfolds independently within each of the lin- 
eages. We calibrated three OU models simulating evolu- 
tion with one (OUl), three (OU3), or five (OU5) 
phenotypic optima (Butler and King, 2004) (see Fig. 5). 
Models were fitted over the complete set of species, as 
well as on subsets of species including only herbaceous 
(forbs and graminoids) or woody (shrubs and trees) 
growth forms (see SOM for details on growth forms and 
model fitting). 

Finally, we also considered a nuU model of phylogeneti- 
cally independent evolution (PI) to test for the existence 
of phylogenetic structure in the trait data. The PI model 
ignores phylogenetic relatedness across species as if they 
were placed at equal distance on a star phylogeny. It 
serves as a null hypothesis which supposes the absence of 
phylogenetic covariance in the trait distribution. 

Treewise phylogenetic signal 

We calculated three statistics to quantify the magnitude 
and test the significance of the overall phylogenetic signal 
in the phylogeny of vascular plants (see SOM for details): 
Pagel's ?. (Pagel 1999), Blomberg's ^T-statistic (Blomberg 
et al. 2003), and phylogenetic autocorrelation measured 
by Moran's index (Gittleman and Kot, 1990). These sta- 
tistics provide information on the likelihood of the phylo- 
genetic pattern with respect to models of trait evolution. 
In Pagel's approach, a value of A = 0 indicates evolution 
independent of phylogeny (PI), whQe a value of 1=1 
indicates that the phylogenetic pattern conforms to 
Brownian motion on the given phylogeny (Freckleton 
et al. 2002). The ^T-statistic measures the degree to which 
a phylogenetic tree correctly describes the covariance 
structure observed in the data compared to a BM model 
along the candidate tree: K > 1 indicates higher conserva- 
tism than in the BM case. Finally, Moran's index, I, quan- 
tifies the similarity in LMA among species with respect to 
their phylogenetic distance (Gittleman and Kot, 1990; 
Diniz-FiUio 2001). It is a measure of the covariation of 
trait values across species weighted by the phylogenetic 
distance between them: / varies from 1, indicating strong 
positive association (similarity across close relatives) to —1, 
indicating strong negative association. It has an expected 
value of where n is the number of tips, under the 
null hypothesis of no correlation between trait values and 



phylogenetic distance. We also calculated / within different 
classes of divergence time (i.e., age of the most recent com- 
mon ancestor) which were chosen to ensure a sufficient 
number of observations within each class. This was used to 
build a resulting phylogenetic correlogram, which repre- 
sents how trait similarity among species varies with the 
time since their divergence. In theory, this temporal pattern 
allows one to discriminate the BM from the OU models 
(Hansen and Martins, 1996; Diniz-Filho 2001): a linear 
decrease in similarity with time conforms to a model of 
evolution by Brownian motion, whereas an exponential 
decrease indicates evolution by drift and stabilizing selec- 
tion as in the OU models (Hansen and Martins, 1996). 
Thus, we adjusted linear and exponential fits to the 
observed empirical patterns using nonlinear least squares 
and compared the fits using a likelihood ratio test. Similar 
analyses were also performed within subsets corresponding 
to the different growth forms. 

Phylogenetic signal at clade level 

We analyzed the phylogenetic signal at clade level using 
the analysis of traits procedure (AOT Ackerly, and Kem- 
bel, 2007). This analysis can detect functionally diversify- 
ing and conservative ancestral splits in trait values across 
daughter clades (Ackerly and Nyffeler, 2004; Moles et al. 
2005). At each divergence in the phylogeny (i.e., node of 
the tree), ancestral mean trait values are calculated 
according to the procedure described in Felsenstein 
(1985). We then calculated the divergence width statistic, 
DW, which measures the degree of divergence (or trait 
radiation) between the child clades of each node. An 
ancestral node was considered as diversifying when the 
corresponding DW value was higher than expected at ran- 
dom, and as conservative when DW was lower than 
expected at random. Significant divergence widths were 
detected by testing the statistic against the null hypothesis 
PI using permutations (n = W^) and correcting for multi- 
ple comparisons (Benjamini and Hochberg, 1995). Signifi- 
cantly diversifying nodes, that is, with higher DW than 
expected under the null hypothesis (PI), indicate past 
diversifying events (trait radiation) followed by subse- 
quent conservatism within daughter clades. On the other 
hand, significantly conservative nodes correspond to 
divergences of low amplitude between daughter clades 
thus displaying similar trait values and phylogenetic con- 
servatism thereafter. 

Results 

Across the complete dataset, the mean LMA value was 
found to be 75.6 g.m^^, with values ranging from 
8.8 g.m^^ in Impatiens parviflora {Balsaminaceae) to 
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1479.1 g.m^ in Hakea leucoptera (Proteaceae). LMA 
ranges largely overlapped across the major clades in 
Tracheophyta, but on average, Lycopodiophyta and 
Monilophyta (ferns) exhibited the lowest LMA, while 
Gymnospermae showed the highest LMA compared to the 
other major clades (Fig. 1; Table SI). Within the Angio- 
spermae, Monocotyledoneae tended to have lower LMA 
values than Eudicotyledoneae (Fig. 1). Finally, within Eudi- 
cotyledoneae (80% of sampled species; Table SI), basal 
clades showed higher LMA values than more recently 
diverged clades (Asteridae and Rosidae; Figs 1 and 2). 

At the family level, within angiosperms, the proportion 
of woody species in a family explained 24% of the varia- 
tion in family mean LMA values (Fig. 3A). This was due 
to the tendency for woody species to have higher LMA 
than herbaceous species (Fig. 81). Nevertheless, the mean 
LMA of herbaceous and woody species belonging to the 
same family was found to be positively related (Fig. 3B). 
This is a first indication of the existence of phylogenetic 
conservatism in LMA values within families which seems 
to interact with the overall effect of growth forms on LMA. 

Overall, common statistics characterizing the phyloge- 
netic signal in LMA indicated low but significant similar- 
ity between close relatives, hence low conservatism, as 



Lycopod. (7) 
IVIonilophyta (74) 
gymno. (81) 
b.Angio. (5) 
b.lVlonocot. (211) 
Comm. (544) 
Magnoliidae (134) 
b.Eudicot. (733) 
Rosidae (1959) 
Asteridae (1653) 



100 
LMA (g m-2) 

Figure 1. Leaf mass per area (LMA) within clades. Distribution of 
LMA (g.m^^) witli respect to the major monophyletic clades of 
vacular plants (Tracheophyta) in increasing order of divergence time 
(letter "b" refers to basal clades, numbers in bracl<ets indicate the 
number of species in the clades, n). Boxes and vertical lines indicate 
the interquartile range and the median in each category. Box height is 
proportional to \/n (see Table SI). Dotted lines represent the 
smoothed distribution within each category, and crosses indicate 
outliers. 




well as little departure from the expected pattern of evo- 
lution following BM: Blomberg's JC-statistics (Blomberg 
et al. 2003) was 0.093 (= ^), Pagel's / (Pagel 1999) was 
estimated to 0.946, and Moran's phylogenetic index (Git- 
tleman and Kot, 1990) was 0.034 (see Material and Meth- 
ods for interpretations). 

In order to identify specific clades showing significant 
evidence of conservatism or diversification, we estimated 
the divergence width, DW, or magnitude of divergence in 
LMA values across daughter clades at each evolutionary 
node (Ackerly and Nyffeler, 2004; Moles et al. 2005). 
When the 1675 internal nodes of the phylogeny were con- 
sidered, DW significantly differed from random in just 62 
cases (3.7%). Of these 62 nodes, 12 corresponded to 
diversifying nodes with trait radiation (Table 1, Fig. 2), 
and 50 to conservative nodes with trait clustering across 
daughter clades. Conservative nodes mostly occurred 
among terminal nodes (genera. Table S4), indicating con- 
servatism in LMA values among congeneric species, and 
at the basis of four larger clades including the Lecythida- 
ceae and Sapotaceae families (Table 1). Cladogenesis at 
the base of these four clades led to sister clades more 
homogeneous than expected at random with regard to 
LMA. Compared to splits showing conservatism in LMA, 
diversifying nodes were typically older (70.2 Myr old vs. 
17.2 Myr old on average), and thus account for much of 
the evolution of the basic range of leaf functioning. The 
most extreme diversifying events in LMA were ancient 
divergences followed by subsequent conservatism within 
descending clades (Table 1, Fig. 2). Diversification in 
LMA thus seems to have occurred early during the course 
of vascular plant evolution. We investigated whether this 
pattern was driven by the contrast between the large 
majority of angiosperms in the dataset and the early- 
diverging lycophyte, fern, and gymnosperm clades by run- 
ning the analysis on the dataset restricted to angiosperms. 
When early- diverging clades were excluded, the same 
early evolutionary splits were identified as diversifying 
within the angiosperm clade (Table 1). 

Phylogenetic correlograms further revealed how simi- 
larity in LMA among species rapidly decreases with diver- 
gence time (Fig. 4). On average, the older the most 
common ancestor, the lower the correlation in LMA tip 
values. Similar values for recent divergences (low diver- 
gence time) showed that, on average, any two close 
relatives among woody species were closer functionally 
than any two close relatives among herbaceous species. 
More generally, woody species showed a greater overall 
similarity through time and hence a higher degree of con- 
servatism in LMA than herbaceous species (Fig. 4). Angio- 
spermae, representing 97% of sampled woody species 
(Table SI), were mostly responsible for this pattern. 
Moreover, functional similarity in LMA decreased differ- 
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Selaginella (1) 
Lycopodiaceae (6,3/4) 
Ophioglossaceae (2,2/3) 
Equisetum (7) 
Osmunda (3) 
Dicranopteris (1) 
Pilularia (1) 
Salviniaceae (2,2/2) 
Cyathea (2) 
Dicksoniaceae (4,3/7) 
Pteridaceae (4,3/40) 
Dennstaedtiaceae (3,2/18) 
Dryopteridaceae (14, 3/47, 1c) 
Polypodiaceae (9,4/79) 
Oleandraceae (2,2/4) 
Blechnum (4) 
Woodsiaceae (9,6/20) 
Asplenium (7) 
Cycas (1) 
Macrozamia (2) 
Ginkgo (1) 
Pinaceae (36,7/13) 
Ephedra (2) 
Araucaria (1) 
Podocarpaceae (19,8/18) 
Taxus (2) 

Cupressaceae (17,6/29) 
Nymphaeaceae (2,2/6) 
lllicium (1) 
Ceratophyllum (1) 
Acorus (1 ) 
Araceae (12,9/109) 
Alismataceae (6,4/12) 
Butomus (1) 

Hydrocharitaceae (7,5/20) 
Scheuchzeria (1) 
Triglochin (2) 

Potamogetonaceae (19,4/7) 
Zostera (1) 
Freycinetia (1) 
Carludovica (1) 
Dioscoreaceae (2,2/8) 
Melanthiaceae (9,5/24) 
Colchicum (1) 
Bomarea (1) 
Liliaceae (9,5/11) 
Smilax (4) 
Philesiaceae (2,2/7) 
Ripogonum (1) 
Hypoxis (1 ) 

Orchidaceae (42,21/1337, 1c) 
Iridaceae (15,5/89) 
Xanthorrhoeaceae (1 1 ,4/1 ) 
Asparagaceae (35,18/1) 
Alliaceae (25,6/31) 
Arecaceae (16,13/204) 
Anigczanthos (1) 
Commelina (2) 
Pontederiaceae (4,3/9) 
Heliconia (1) 
Costus (1 ) 

Zingiberaceae (2,2/47) 
Marantaceae (3,3/32) 
Bromeliaceae (4,2/51) 
Typhaceae (5,2/2) 
Cyperaceae (160, 18/1 02, 1c) 
Juncaceae (42,2/7) 
Flagellaria (1) 
Restionaceae (3,3/40) 
Poaceae (299,104/707,2c) 
Hedyosmum (2) 
Winteraceae (6,3/4) 
Piperaceae (11,2/10) 
AristDlDchiaceae (3,2/8) 
Hernandia (1) 
Monimiaceae (10,8/39) 
Lauraceae (64,19/49) 
Myristicaceae (8,6/18) 
Magnoliaceae (7,3/7) 
Eupomatia (1) 
Annonaceae (23,17/125) 
Eudicots 



Saxifragalesl 



Crossosomatale: 



Sapindalk 




Malpighialei 





(C) 

Papaveraceae (18,8/41) 
Euptelea (1) 
Decaisnea (1) 
Menispermaceae (3,3/73) 
Ranunculaceae (87, 15/52, 1c) 
Berberidaceae (10,4/18) 
Platanus (2) 

Proteaceae (193, 29/69, 3c) 
Meliosma (6) 
Tetracentron (1) 
Buxus (2) 
Gunnera (1) 

Hamamelidaceae (7,4/29) 
Daphniphyllum (1) 
Cercidiphyllum (1) 
Ribes (12) 

Saxifragaceae (17,3/32) 
Paeonia (2) 
Haloragaceae (4,3/9) 
Crassulaceae (22,6/39) 
Vitaceae (10,5/14) 
Geraniaceae (21 ,2/1 1,1c) 
Stachyurus (1) 
Staphyleaceae (3,2/5) 
Combretaceae (18,6/20) 
Lythraceae (10,7/30) 
Onagraceae (32,6/18) 
Myrtaceae (146,36/127,1c) 
Vcchysiaceae (5,3/7) 
Olinia (1) 

Melastomataceae (21 , 1 0/1 94) 
Tropaeolum (1 ) 
Reseda (5) 

Brassicaceae (125,47/426, 1c) 
Malvaceae (72,38/273) 
CochlDspermum (2) 
Thymelaeaceae (13,6/58) 
Cistaceae (23,4/8) 
Dipterocarpaceae (9,4/16) 
Sapindaceae (55, 18/1 38, 1c) 
Anacardiaceae (30,17/68) 
Burseraceae (8,4/18) 
Rutaceae (40,28/162) 
Meliaceae (22,14/51) 
Simaroubaceae (6,5/24) 
Zygophyllaceae (16,8/28) 
Fabaceae (363,112/677,2c,2d) 
Polygalaceae (10,2/19) 
Rosaceae (219,37/1 04,4c) 
Elaeagnaceae (5,2/3) 
Rhamnaceae (35,12/51) 
Ulmaceae (31,7/15) 
Cannabaceae (2,2/2) 
Urticaceae (26,12/49) 
Moraceae (39, 13/37, 1c) 
Cucurbitaceae (6,4/121) 
Coriaria (2) 
Fagaceae (91,8/8,1c) 
Juglandaceae (15,4/8) 
Myricaceae (7,2/3) 
Ticodendron (1 ) 
Casuarinaceae (7,2/4) 
Betulaceae (42, 6/6, 2c) 
Parnassia (2) 
Celastraceae (25,9/86) 
Cunoniaceae (13,8/23) 
Elaeocarpaceae (10,5/9) 
Connarus (2) 
Oxalis (6) 
Cordobia (1) 
Elatine (2) 

Euphorbiaceae (115,49/333, 1c) 
Linaceae (8,2/14) 
Ochnaceae (3,2/26) 
Irvingiaceae (2,2/3) 
Chrysobalanaceae (7,3/17) 
Ixonanthes (1) 
Passiflora (1 ) 
Rhizophoraceae (10,5/14) 
Clusiaceae (34,12/47) 
Caryocar (2) 
Violaceae (27,5/20) 
Goupia (1) 
Microdesmis (1 ) 
Humiria (1) 
Lozania (1 ) 

Salicaceae (88,14/81 ,2c) 
Clade C 



Caryophyllaled 




Aqui foliales . 



Dipsacalesj^ 



— O^steralBi 



Gentranales 



Lamiat^s 



Aextoxicon (1 ) 
Dilleniacsae (16,4/11) 
Loranthaceae (7,4/77) 
Santalaceae (10,6/37) 
Olacaceae (5,4/27) 
Drosera (3) 
Myricaria (2) 

Polygonaceae (58,1 1/49,1c) 
Plumbaginaceae (5,3/25) 
Simmondsia (1) 
Amaranthaceae (99,24/183) 
Caryophyllaceae (98,20/93,4c) 
Corrigiola (1) 
Cactaceae (8,8/101) 
Portulacaceae (9,6/25) 
Nyctaginaceae (7,6/38) 
Phytolacca (3) 
Alzoaceae (9,8/128) 
Cornaceae (17,5/8) 
Hydrangeaceae (6,4/16) 
Impatiens (6) 
Polemonlum (1) 
Sapotaceae (25,16/59) 
Maesa (1) 
Jacquinia (1) 
Primulaceae (26,6/22) 
Myrslnaceae (27, 10/42, 1c) 
Lecythldaceae (13, 7/25, 1c) 
Ebenaceae (18, 2/3, 1c) 
Theaceae (14,9/25) 
Symplocos (9) 
DIapensiaceae (2,2/7) 
Styracaceae (6,2/1 1 ) 
Sarracenia (1) 
Clethra (8) 
Cyrilla (1) 

Ericaceae (90,30/1 53,2c) 
Aqulfoliaceae (17,2/2) 
Helwlngia (1) 

Caprlfoliaceae (47, 11/11, 1c) 
Adoxaceae (23,3/5) 
Escallonlaceae (7,7/20) 
Grisellnia (2) 
Arallaceae (22,14/48) 
Aplaceae (99,55/439) 
Plttosporaceae (7,3/9) 
Stylldlum (1) 

Campanulaceae (46,13/85) 
Menyanthaceae (2,2/5) 
GoDdenlaceae (11,4/12) 
Acicarpha (1) 

Asteraceae (479,147/1511,6c) 
Aucuba (1) 
Icaclnaceae (4,4/53) 
Boraglnaceae (55,30/151) 
Convolvulaceae (13,7/56) 
Solanaceae (33,11/97) 
Rublaceae (97,41/620) 
Apocynaceae (43,25/496) 
Gentlanaceae (33,8/77) 
Loganlaceae (8,4/20) 
Oleaceae (36,13/26) 
Plantaglnaceae (85,20/23) 
Scrophularlaceae (30,7/263, 1c) 
Lentibularlaceae (4,2/3) 
LIndernIa (1) 
MImulus (3) 

Orobanchaceae (33, 8/65, 2c) 
Paulownia (1) 
Acanthaceae (1 1 ,8/235) 
BIgnonlaceae (1 6,14/1 1 1 ) 
Verbenaceae (18,11/85) 
Lamiaceae (91 ,35/213, 1c) 
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Figure 2. Phylogenetic tree of the sampled extant species. The tips of the tree correspond to botanical families, except for families with one 
genus only, in which case the genus is indicated. Symbol size represents the clade mean LMA in five classes of increasing values: < BO, [50-75], 
[75-100], [100-150], and > 150 g.m"^. Red symbols indicate significant conservative evolutionary splits (with divergence width, DW, 
significantly lower than expected at random), green symbols represent significant diversifying divergences (DW, significantly higher than expected 
at random). Numbers in parentheses are: number of sampled species, number of sampled genera / number of genera in the family, according to 
Kew's classification. Where appropriate, the number of significant divergences within families is indicated in parentheses as "c" for conservative 
and "d" for diversifying. Tracheo.: vascular plants (Tracheophytes), Euphyllo.: megaphyll plants, Spermato.: seed plants (Spermatophytes). Branch 
lengths are indicative. 




0 20 40 60 80 100 50 100 200 500 

% woody species LMA woody (g m-^) 



Figure 3. Phylogenetic and growth form signal in LMA across and within angiosperms families: (A) mean LMA, LMA, versus proportion of woody 
species (shrubs and trees; best linear fit: r' = 0.24, P < 10"*); (B) Mean LMA of herbaceous species (forbs and graminoids) versus mean LMA of 
woody species within families (dashes: 1 to 1 line; best linear fit: r^ = 0.25, P < 10"^). Symbols and colors indicate major clades of angiosperms: 
O Eudicotyledoneae (orange: basal families, red: Asteridae, dark red: Rosidae); a Monocotyledoneae (green: basal families; dark green: 
Commelinidae); ^ Magnoliidae. Only families with over 10 sampled species were considered. 



Table 1. Significant evolutionary splits below genus level (see Table S4 for results concerning genera). Clades are named after taxonomic names, 
unless undefined. Sister clades, separated by "/", descend from the focal nodes (not detailed for the polytomies at the base of the Lecythidaceae 
and Sapotaceae families); Xantho. indicates the Xanthorreaceae family. Age: estimated age of the evolutionary split in Myr. tt tips: number of 
descending tips (size of the clade). LMA: respective mean LMA value within sister clades, in g.m"^, except for Lecythidaceae and Sapotaceae 
(mean LMA within parent family). Woody (%): respective proportion of woody species (shrubs/trees) within both sister clades. DW: divergence 
width of the split. Diversifying and conservative splits showed significantly higher and lower DW compared to random, respectively. 



Name 


Sister clades 


Age (Myr) 


tt tips 


LMA (g.m-^) 


Woody % 


DW 


Diversifying splits 














Eudicotyledoneae 


Ranunculales 1 rest of Eudicotyledoneae 


144 


4345 


51/78 


11/55 


0.598 




Proteales 1 sister clade 


142 


4225 


301/73 


1 00/53 


0.619 


Asparagales 


Alliaceae + Asparagaceae + Xantho. 1 Orchldaceae + Hypoxidaceae 


108 


129 


94/43 


6/0 


0.682 


Ericales 


Impatlens 1 rest of Ericales 


91 


249 


17/96 


0/81 


1.316 




Alliaceae + Asparagaceae 1 Xantho. 


74 


71 


72/287 


3/27 


0.931 




Cyperaceae 1 Juncaceae 


61 


202 


66/83 


0/0 


0.528 




Scrophulariaceae 1 9 families of Lamiales 


55 


208 


101/56 


53/29 


0.662 




Mirbelieae tribe / sister clade in Fabaceae 


53 


192 


1 45/52 


1 00/1 7 


0.639 


Brassicaceae 


Capparoideae subfamily / rest of Brassicaceae 


43 


125 


161/44 


100/2 


1.081 




Betulaceae 1 Casuarinaceae 


30 


49 


65/445 


1 00/1 00 


1.442 




Portulacaceae 1 Cactaceae 


22 


17 


66/523 


0/13 


1.467 




Acacia 1 Ingeae tribe in Fabaceae 


20 


63 


1 95/83 


100/100 


0.659 


Conservative splits 
















Ebenaceae + Lecythidaceae 1 sister clade 


69 


162 


1 02/1 1 3 


97/96 


<io-^ 


Sapotaceae 


1 6 sister clades 


45 


25 


124 




0.195 


Lecythidaceae 


7 sister clades 


35 


13 


99 


100 


0.099 




Sagina 1 Bufonia tenuifolia in Car/o. 


16 


7 


63/63 


0/0 


<10"^ 
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Figure 4. Similarity in LIVIA as a function of divergence time. Moran's index (Im) in classes of divergence time (IVIyr) for (A) tiie complete species 
set (n = 5401), and (B) tiie iierbaceous (□, n = 2417) and woody (A, n = 2564) growtii forms. Filled (resp. open) symbols indicate (non-)significant 
correlations (tzievel: 0.05). Error bars were smaller than symbol size and tinerefore not represented. Dotted lines represent linear fits (Brownian 
IVIotion model, BIVI); dashed lines represent exponential fits (Ornstein-Uhlenbeck model, OU). Arrows indicate major evolutionary nodes; Tracheo.: 
vascular plants, Euphyllo.: megaphyll plants (ferns + seed plants), Sperma.: seed plants, gymno.; Gymnospermae, Anglo.; Angiospermae, IVIagno.; 
Magnoliidae. 



ently in woody species (shrubs and trees) and in herba- 
ceous species (forbs and grasses; Fig. 4B), suggesting that 
resource-use strategies have evolved conjointly with 
growth form in vascular plants. 

We further analyzed the evolution of LMA by means of 
models of continuous trait evolution. The models sup- 
ported the hypothesis of stabilizing selection around some 
optimal trait values against evolution of LMA by drift 
only: the OU performed better than the simple BM mod- 
els according to BIC ranking (Table S3). This conclusion, 
drawn for the complete phylogeny, was also supported 
for woody and herbaceous growth forms considered sepa- 
rately. The better fit of OU models was consistent with 
the linear trends of decreasing similarity with divergence 
time observed in the empirical correlograms (Fig. 4), as 
opposed to an expected exponential decrease in the case 
of BM models (Diniz-FUho 2001). Second, models incor- 



porating different selection regimes (Fig. 5), that is, 
stabilizing selection around different optimal trait values 
(Butler and King, 2004) performed better than single- 
optimum models (Table S3), supporting the existence of 
shifts in LMA selection regimes between Angiospermae 
versus Gymnospermae, and between Monocotyledonae ver- 
sus Eudicotyledonae. Finally, such shifts in selective 
regimes were also observed within each growth form ana- 
lyzed separately: between the Angiospermae and Gymno- 
spermae clades for woody species, and between the 
Monocotyledonae and Eudicotyledonae clades for herba- 
ceous species (Fig. 5, Table S3). The best identified mod- 
els explained, respectively, 26.3%, 25.9%, and 42.3% of 
the null deviance for the complete, woody, and herba- 
ceous datasets. 

Orstein-Ulhenbeck models allow us to estimate the rate 
of adaptation to phenotypic optima in LMA (a parameter 




BM (■) 

or 

0U1 (•) 




0U3(A) 




OU5(T) 



Lycopodiophyta (7) 
Monilophyta (74) 
Cycadophyta (3) 

Ginkgo (1) 

Coniferae (77) 

Nymphaeaceae (2) 
Schisandraceae (1 ) 
C9ratophyllac9a9 (1 ) 

Monocotyledonae (755) 

Chloranthaceae (2) 
Magnoliidae (134) 
Eudicotyledonae (4345) 



Figure 5. Alternative models of LIVIA evolution based on Brownian motion (BM, ■) and Ornstein-Uhlenbeck process with one (0U1, •), three 
(0U3, A ). 01" five (0U5, ▼ ) phenotypic optima indicated by different colors. Only major and basal clades of Tracheophyta are represented 
(arbitrary branch lengths; numbers indicate clade size). IVIodels were fitted on the complete phylogeny and on the herbaceous and woody species 
groups separately by pruning the phylogeny accordingly (see Fig. S6). 
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on Fig. 6A, Hansen 1997). The quantity t = measures 
the time taken by a phenotype evolving under a new 
regime to move halfway from its ancestral state (Hansen 
1997). We found that this time was 15.5 Myr for herbs 
and 24.5 Myr for woody species, that is, a longer time 
scale along woody lineages for new selective conditions to 
be more influential than constraints from the ancestral 
state. 

Discussion 

The evolutionary history of a key plant functional trait, 
leaf mass per area (LMA), was investigated across a large 
set of vascular plant species. Using comparative phyloge- 
netic methods and a global trait database of over 5000 
species, we detected that leaf strategies were overall 
weakly, albeit significantly, conserved during the diversifi- 
cation of vascular plants. This low conservatism appeared 
to be associated with multiple regimes of stabilizing selec- 
tion, that is differing adaptive optima across the major 
clades. Importantly, we highlight the strong interaction 
between the evolution of growth forms and of leaf strate- 
gies, with woody species displaying slower leaf diversifica- 
tion and higher conservatism in LMA than herbaceous 
plant species. 

Leaf mass per area appeared significantly but weakly 
conserved. This low phylogenetic signal captures the fact 
that distantly related species may occupy similar positions 
along the leaf economics spectrum. The large variability 
found within clades here evidenced a pattern of func- 
tional convergence across clades over a broad phyloge- 
netic scale. However, significant levels of trait 
conservatism imply that, to a certain extent, closely 
related species do tend to invest dry matter similarly 
within leaf tissues and hence occupy nearby positions 



along the leaf economics spectrum. In particular, evolu- 
tionary splits identified as conservative nodes were mostly 
found at the level of genera. As might be expected, the 
degree of similarity in resource-use strategies varied with 
the length of the period of common evolution along lin- 
eages (Hansen and Martins, 1996). This pattern indicates 
a temporal dimension to the diversification of resource- 
use strategies among species. When considering large time 
scales, the phylogenetic pattern in LMA displayed early 
diversification associated with the divergences between 
extant ferns and lycophytes on the one hand and seed 
plants (Spermatophytae) on the other, and between Angio- 
spermae and Gymnospermae (Fig. 4). Within the Angio- 
spermae clade, diversity in LMA reflects the wide adaptive 
radiation of flowering plants into a range of ecological 
strategies, involving diverse growth forms and ecological 
niches both within and across habitats (Ricklefs and 
Renner, 1994; Losos 2008). By contrast, the consistency of 
leaf economics in the other major clades may reflect both 
stabilizing selection (around a given optimum) as well as 
their marginalization to particular ecological situations 
during the adaptive radiation of the Angiospermae. In fos- 
sil floras, basal plant clades display a higher diversity of 
growth forms than is currently observed among extant 
species: treelike plants existed within Lycopodiophyta, as 
well as ruderal herbaceous taxa within gymnosperms 
(Roth well et al. 2000). The fossil record thus reveals that 
extinction events also shaped the spectrum of resource- 
use strategies that exist today in vascular plants. 

Models of continuous trait evolution supported the 
hypothesis that stabilizing selection shaped the phyloge- 
netic pattern in LMA, against the hypothesis of evolution 
by drift alone. Given the breadth of sampling in our 
dataset, over a large taxonomic scale and a range of eco- 
logical situations, the stabilizing selection we detected 
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Figure 6. Parameter estimates for models of trait evolution: (A) parameter a estimates tiie rate of adaptation (Hansen, 1997), (B) parameter a 
estimates the magnitude of perturbations not due to selection (Hansen, 1997). Symbols indicate different models: Brownian motion (■) and 
Ornstein-Uinlenbeck models, witli one (•), three ( A ). oi" five ( ▼ ) selective optima (see Fig. 5 for model description). White symbols indicate 
the best models according to BIC ranking (Table S3). Error bars were too small to be represented. 
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likely reflects selection against extreme phenotypes that 
infringe structural and physiological limitations to LMA. 
The phylogenetic patterns evidenced here showed simulta- 
neously conservatism, mostly within genera, and conver- 
gence across large clades. Congruence between radiation 
and restrictive environmental conditions, such as within 
the Proteales, also suggests strong filtering effects that 
may lead to conservatism at larger phylogenetic scales and 
functional distinctiveness (Cornwell et al. 2014). Physical 
conditions in natura set constraints on plant strategies 
that may directly or indirectly affect achievable LMA val- 
ues. A recent study of the evolution of the leaf economics 
spectrum highlighted the nature of evolutionary forces on 
LMA (Donovan et al. 2011). The authors argued that 
selection in general, rather than genetic constraints, direc- 
ted the evolution of LMA. As an example, direct selection 
could influence LMA via biomechanical limits on the 
amount of dry matter needed to support and maintain a 
planar photosynthetic surface. Indirect selection could 
result from selection against unfit LMA values in relation 
to major trade-offs in plant strategies (Donovan et al. 
2011). Furthermore, the models of trait evolution deter- 
mined different optimal phenotypes across major clades 
of vascular plants. These findings indicate that these 
clades have evolved under distinct selective regimes (But- 
ler and King, 2004) that may be interpreted as ecophysio- 
logical constraints that support different LMA optima. 
These constraints can be understood as trade-offs in plant 
traits that may render particular values of LMA unfit 
(Donovan et al. 2011). Changes in functional optima 
along the phylogeny have likely resulted from a combina- 
tion of changing selection pressure effects, trade-offs 
among axes of plant variation imposed by these pressures, 
and historical contingency or ancestry. These variations in 
LMA optima across clades also existed within growth 
forms. Hence, major divergences during plant evolution 
gave birth to clades that in time evolved leaf strategies 
optimized in response to selective regimes that have 
imprinted the evolution of both woody and herbaceous 
lineages within those clades. 

Recent analyses of molecular sequences have shown that 
woody species have evolved more slowly than herbaceous 
species (Smith and Donoghue, 2008) and taller species 
more slowly than shorter species (Lanfear et al. 2013), pre- 
sumably because of their longer generation times. However, 
the issue of whether molecular evolutionary rates coincide 
with phenotypic evolutionary rates is somewhat controver- 
sial (Bromham et al. 2002). Here, we found evidence sup- 
porting different evolutionary rates for LMA across growth 
forms, with higher rates in herbaceous compared to woody 
species. The woody syndrome occurred with higher and 
more consistent functional similarity among lineages that 
shared the habit. We suspect that allometry in carbon allo- 



cation patterns inherent to woody plants has constrained 
the variability of achievable resource-use strategies and led 
to higher conservatism along woody lineages through time. 
Woody plants must allocate carbon extensively to second- 
ary cell wall thickening (including lignification), particu- 
larly in mechanical tissues and secondary xylem tissue. This 
requisite feature of carbon allocation constrains the 
resource-use strategies of woody plants. Vascular plants in 
general exhibit a three-way trade-off between: small leaves 
of high LMA, large leaves of intermediate LMA, and small 
leaves of low LMA (Pierce et al. 2013). This latter combina- 
tion, common among ruderal herbaceous species and aqua- 
tic plants, is rarely exhibited by woody species. 
Furthermore, many of the extremely high LMA gymno- 
sperm taxa exhibit phylogenetic constraint based on the 
structure and functioning of the xylem: tracheary elements 
of most extant gymnosperm families consist of solely trac- 
heids, and not of tracheae that would provide sufficient 
internal translocation to support fast growth rates and 
extensive evapotranspiration from broad leaves. Gymno- 
sperms therefore typically exhibit a suite of sclerophytic 
traits, such as a thick leaf endodermis, alongside resin 
canals to guard against predator and pathogen attack, 
embodied in "needle" leaves of high LMA. These adapta- 
tions require an investment in carbon and mineral 
resources in leaf mass that cannot then be allocated to low 
LMA leaves typically associated with fast growth rates in 
ruderal angiosperms. Thus, the investment of carbon in 
wood and additionally the type of wood produced both 
provide constraints to woody plant resource-use strategies. 

The observed patterns in the evolution of LMA raise 
the issue of selective pressures resulting in different selec- 
tion regimes at large evolutionary scales. Fossil evidence 
suggests that large spatiotemporal scale changes in climate 
and disturbance regime were partially responsible for the 
diversification of vascular plants (Stebbins 1974; Jacobs 
et al. 1999) and this was likely paralleled by diversifica- 
tion in LMA. Recent evidence showed that plant species 
tend to be more woody and taller toward the tropics 
(Moles et al. 2009). The herbaceous syndrome might have 
evolved under more temperate climatic conditions (Rick- 
lefs and Renner, 1994), which suggests that the woody/ 
herbaceous dichotomy represents a major geographical 
and macroevolutionary dichotomy between climate 
extremes (Zanne et al. 2013). Evolutionary trends in LMA 
probably arose in various ways, as this trait integrates sev- 
eral different aspects of leaf morpho-anatomy (Garnier 
and Laurent, 1994). The interaction of LMA evolutionary 
patterns with growth forms also suggests strong linkages 
between leaf and stem traits in the evolution of resource- 
use strategies (Zanne et al. 2013). 

By exploring the evolutionary history of LMA we have 
revealed different evolutionary trajectories involving coor- 
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dination and trade-offs between leaf traits and plant traits 
related to growth form. As a consequence, different opti- 
mal phenotypes of resource-use strategies have evolved in 
vascular plants setting fundamental limits on the structure 
and physiology of organisms and therefore on ecological 
processes within communities. Continuing advances in 
genomics (Leebens-Mack et al. 2006), phylogenetic infor- 
matics (Smith et al. 2011), and global trait databases 
(Kattge and the TRY consortium, 2011) promise to bring 
further insight into the evolutionary patterns of a range 
of plant functional traits (e.g., plant height, wood density, 
roots traits). Confronting these multiple evolutionary pat- 
terns will provide a powerful general narrative concerning 
the mode and speed of evolution of trait variation and 
distinct plant ecological strategies. 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Appendix SI. References. 

Figure SI. Leaf Mass per Area (LMA) within growth 
forms. Distribution of LMA (g.m^^) with respect to the 
different growth forms in order of increasing mean LMA 
values. In brackets: n, number of species per growth form 
(Table SI). 

Figure S2. Detailed distribution of Leaf Mass per Area 
(LMA) within ferns and lycophytes. 



Figure S3. Leaf Mass per Area (LMA) in the subsample 
of (a) herbaceous species for which life history was unam- 
biguous (annual: mean: 45g.m coefficient of variation: 
/0.51, n = 545; and perennial: 65/0.82, n = 921) and of 
(b) woody species for which leaf phenology was unambig- 
uous (deciduous: 79/0.38, n = 522; evergreen: 174/0.76, 
n = 969). 

Figure S4. LMA at the level of famihes: (a) Mean LMA 
(LMA in g.m^^) across extant species within families vs 
estimated age of families (in Myr); (b) Coefficient of vari- 
ation of LMA within families vs estimated age. 
Figure S5. Correlogram of LMA showing Moran's I (IM) 
within taxonomic levels, for the complete sample (black line 
and symbols, n = 5401 species), and for the woody (dark 
gray, n = 2564 species) and herbaceous (hght gray, n = 2417 
species) growth forms. 

Figure S6. Alternative models for LMA evolution in woody 
(top) and herbaceous species (bottom). 
Table SI. Growth forms and clades. 

Table S2. Results from family-level linear models of the 
mean and coefficient of variation (CV) of log-transformed 
mean LMA values, 1, within Eudicotyledonae families. 
Table S3. Summary of evolution models for LMA. 
Table S4. Genera associated with significant conservative 
divergences. 
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